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\q ABSTRACT 

O 
O . 

The recent massive measurements of pulsar frequency second derivatives have 
shown that they are 100-1000 times larger than expected for standard pulsar 
slowdown low. Moreover, the second derivatives as well as braking indices are 
even negative for about half of pulsars. We explain these paradoxical results on 
the basis of the statistical analysis of the rotational parameters u, v and v of 
the subset of 295 pulsars taken mostly from the ATNF database. We have found 
strong correlation of v and z> either for i) > (correlation coefficient r « 0.9) 



> 

in , 

t^- . and i) < (r « 0.85), and of v and v (r « 0.7). We interpret these dependencies 

as evolutionary ones due to z> being nearly proportional to characteristic age 



c3 



r c h. The derived statistical relations as well as "anomalous" values of v are 
well explained in the framework of the simple model of cyclic evolution of the 
rotational frequency of the pulsars. It combines the secular change of u tr (t), 
Vtr(t) and Vtr{t) according to the power law with n »s 5 and harmonic oscillations 
of 100-1000 years period with an amplitude from 10~ 3 Hz for young pulsars to 
10~ 10 Hz for elder ones. The physical nature of these cyclic variations of the 
rotational frequency may be similar to the well-known red timing noise, however, 
with much larger characteristic time scale. 

t 

Subject headings: methods: data analysis — methods: statistical — pulsars: 
general 



1. Introduction 

The slowdown of radio pulsars is due to the electromagnetic emission energy losses. 
According to this "classical" approach, their rotational frequencies v evolve following the 
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slowdown law i> = —Ku n , where K is the positive constant that depends on the magnetic 
dipole moment and the moment of inertia of the neutron star, n is the braking index. The 
latter can be determined observationally from measurements of u, v and v as n = vv j v 2 . For 
pure dipole structure of pulsar magnetosphere n — 3, the pulsar wind decreases this value 
down to n — 1, while for multipole magnetic field it is up to n = 5 (Manchester & Taylor 1977; 
Manchester, Durdin & Newton 1985; Blandford & Romani 1988). For such law of pulsar 
evolution the second derivatives of their rotational frequencies (and therefore their braking 
indices) can't be measured for objects with the characteristic age r ch = — z//2z> > 10 6 years 
as their v < 1CT 30 s~ 3 , that is much lower than the current detection limit (Livingstone et al. 
2005). At the same time the massive measurements of pulsar frequency second derivatives v 
have shown that their values are much larger than expected for standard spin-down law and 
are even negative for about half of all pulsars. The corresponding braking indices are from 
-10 6 to 10 6 (D'Alessandro et al. 1993; Chukwude 2003; Hobbs et al. 2004). 

These strange results are usually interpreted as an influence of a red timing noise, which 
manifests itself as irregular variations of pulse arrival times relative to the secular frequency 
evolution (Groth 1975; Cordes & Helfand 1980; Cordes & Downs 1985), described through 
Taylor series expansion as 

v(t) = z/ + z> (t - t ) + ^i) (t - t ) 2 + ... (1) 

where to, u , i> and i) are initial time moment, frequency and its derivatives at this moment. 
Historically, the secular evolution of pulsar rotation rate has been considered as the linear 
trend of the frequency behaviour, while quadratic term of the series (1) (with cubic behaviour 
in pulse arrival times), measured as £>, has been treated as the signature of non-monotonic 
(on time scale of the observations) processes - red noise, post-glitch recovery, pulsar orbital 
motion and free precession (Groth 1975; Cordes & Helfand 1980; Shemar & Lyne 1996; Stairs, 
Lyne & Shemar 2000). However, recently the increase of the sensitivity of the receivers and 
time spans of the continuous timing (from several to thirty years) made it possible to estimate 
the higher-order coefficients of the series (1) (Camilo et al. 2000; Scott, Finger & Wilson 2003; 
Hobbs et al. 2004; Livingstone et al. 2005). The results of these studies show that i), while 
changing on the time scales shorter than observation time spans, on average characterizes 
the pulsar secular evolution rather than noise processes. The latter, causing the variations 
mentioned earlier, worsens the accuracy of the second derivative determination, while can't 
change its mean value hundreds times (Camilo et al. 2000; Hobbs et al. 2004; Livingstone et 
al. 2005). 

In this Letter we report the detection of correlation of is and i> either for i) > and i) < 
and correlation of v and v. We show that these statistical relations as well as "anomalous" 
values of v can be explained by the combination of secular decrease of pulsar frequencies 
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Fig. 1. — The i) — z> diagram. The figure shows the pulsars from the work (Hobbs et al. 2004) 
as circles, and the objects measured by other groups as squares. Empty symbols represent 
the pulsars associated with the supernova remnants, and therefore - relatively young ones. 
Analytical fits for both positive and negative branches are shown as solid lines. 



with n 5 and their cyclic variations on the time scale of several hundreds of years. 



2. Statistical analysis of the ensemble of pulsars 

Our statistical analysis is based on the assumption that numerous measurements of the 
pulsar frequency second derivatives reflect their secular evolution on the time scale larger 
than the duration of observations, and uses the parameters of nearly 300 pulsars. 

From 389 objects of ATNF catalogue (Manchester et al. 2005) with known i) we have 
compiled a list of "ordinary" radio pulsars with P > 20 ms and P > 10~ 17 s s _1 , excluding 
recycled, anomalous and binary pulsars, and with relative accuracy of second derivative 



-4- 




measurements better than 75%. It has been appended with 26 pulsars from other sources 
(D'Alessandro et al. 1993; Chukwude 2003). The parameters of all pulsars have been plotted 
on the i) — v diagram (Fig. 1). Here the circles are the objects with second derivatives 
measured at the Jodrell Bank Observatory (Hobbs et al. 2004), squares - the results of other 
groups (D'Alessandro et al. 1993; Chukwude 2003), and empty symbols represent young 
pulsars associated with supernova remnants. 

The main result of the statistical analysis of these data is a significant correlation of v 
and z>, either for 168 objects with i) > (correlation coefficient r pa 0.90) and for 127 objects 
with v < (r pa 0.85). Either groups follow the nearly linear laws, however they are not 
exactly symmetric with respect to v — 0. Also, the young pulsars with z> < — 10~ n s~ 2 and 
i> < are absent. 

Figure 2 shows the dependence of v on the characteristic age r c h- Of course, due to 
practical proportionality of these parameters there is a strong dependence of v on r c ^ either 
for pulsars with v > (positive branch, r pa 0.85) and for v < (negative branch, r pa 0.75) 
(Fig. 3a). Also, we have found the significant correlation of n and r ch (Fig. 3b). Note some 
important points of these results. 



• Measurements of different observational groups show the same laws, so the correlations 
discovered are not related to the peculiarities of the receiving devices of a single radio 
telescope. 

• Young pulsars confidently associated with supernova remnants are systematically shifted 
to the left in Figure 1 and Figure 2 (empty symbols). The order of their physical ages 
roughly corresponds to that of characteristic ages. It means that any dependence on 



- 5 - 



v or r c h may be interpreted as the dependence on time. 

• Separate analysis of the whole set of parameters of pulsars with v > and i) < (galac- 
tic positions, fluxes, luminosities, magnetic fields, dispersion and rotation measures) 
shows the statistical identity of these sub-groups. 

So, the 'v—v diagram (Fig. 1) may be interpreted as an evolutionary one. In other words, 
each pulsar during its evolution moves along the branches of this diagram while increasing 
the value of its v (which corresponds to the increase of its characteristic age). However, 
there is an obvious contradiction, as for i) < (negative branch), z>, being negative, may 
only decrease with time, so the motion along the negative branch may be only backward! 
Certainly, this contradiction is easily solved by assuming non-monotonic behaviour of v(t), 
which may have some cyclic (quasi-periodic) component along with secular one. Figure 4 
shows the example of such cyclic behaviour of a pulsar. The characteristic time scale T 
of such cycles must be much shorter than the pulsar life time and at the same time much 
larger than the time scale of the observations. Each pulsar during its evolution repeatedly 
changes the sign of i) and spends roughly half of its life time on each branch. The asymmetry 
of the branches reflects the positive secular behaviour of v tr (t), and so, secular increase of 
v{t). Systematic decrease of branches separation reflects the decrease of the oscillations 
amplitude and/or the increase of its characteristic time scale. Any non-monotonic variations 
of v(t), like post-glitch recovery, standard red noise or orbital motion, will manifest itself 
in a similar way on the v — v diagram and lead to extremely high values of i) (Groth 1975; 
Cordes & Helfand 1980; Shemar & Lyne 1996; Stairs, Lyne & Shemar 2000). However, their 
characteristic time scales vary from weeks to years, and they are detected immediately. In 
our case we have processes on time scales of hundreds of years, and their study is possible 
only by means of statistical methods, assuming that the ensemble of pulsars of various age 
models the behaviour of a single pulsar during its evolution. 

Note that the correlations between \P\ and P, similar to these discovered here (but 
without distinguishing between objects with i) > and i) < 0), were discussed earlier (Cordes 
& Downs 1985; Arzoumanian et al. 1994; Lyne 1999) as an evidence of higher timing noise 
of young pulsars. 

3. The simplest model of pulsar cyclic evolution 

Cyclic variations of the pulsar rotational frequency may be complicated - quasi-periodic, 
multi-harmonic, stochastic. As a rough initial approximation, we use a simple model of 
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harmonic oscillations superimposed on secular evolution. In this case, 

u{t) = u tr {t) + A(t) sin [Q(t)t + O ] , (2) 

where A(t) is the amplitude of periodic frequency variations, Q(t) = 2ix/T{t) - its frequency, 
and v tr {t) is the secular evolution term. This leads to 

v{t) = v tr {t) + A{t)VL{t) cos [n(t)t + O ] , (3) 

j){t) = is tr (t) - A(t)tt 2 (t) sin [Sl(t)t + <j> Q ] , (4) 

where we neglected the terms with A(t) and Q(t) derivatives due to its relative smallness. 
From (3) and (4) it is easy to get the relation between i)(t) and v(t) in the cycloidal form 
[i){t) - i>tr(t)} 2 = Vl 4 (t)A 4 (t) - Vl 2 {t)[v{t) - u t r(t)] 2 , which is plotted in Figure 4. 

By marking with "+" and "-" subscripts the values related to the positive and negative 
branches on the v — v diagram (Fig. 1), correspondingly, and pointing to the amplitude 
values of v and v (i.e. assuming sin [Q(t)t + O ] = Tl), w e get 

u ± (t) = u tr (t) t A(t) (5) 

U ± (t) = Vtr{t) (6) 

i) ± (t) = i) tr (t)±A(t)n 2 (t). (7) 

So, due to (6), the branches of v — u diagram describe the dependence of pulsar frequency 
second derivative v on the secular behaviour of its first derivative v tr , which corresponds to 
the dependence on time. 

Relations (2) - (7) describe the apparent evolution of some "average" pulsar. Of course, 
in the ensemble of all pulsars the scatter of points on the v — v diagram reaches 4 orders of 
magnitude due to distribution of the neutron star parameters and initial conditions. Due to 
this fact it is impossible to determine Utr{t) directly as the half-sum of v+(t) and i>-{t) from 
(7) and Fig. 1, as it requires the subtraction of two close values with significant errors. 

First of all we compute the behaviour of vtr{t) (or, v tr [ytr)) by plotting the studied 
pulsar group on the v — v diagram (Fig. 5). The objects with i) > and v < are marked 
as filled and empty circles, correspondingly. It is easily seen that the behaviour of these two 
sub-groups is the same, i.e. A(t) according to (5) is significantly smaller than the intrinsic 
scatter of v{y). However, a strong correlation between u tr and u tr (r 0.7) is seen, and 

hr = -Cvl, (8) 

where C = icr 15 - 26±L38 and n = 5.13 ± 0.34. So, the secular evolution of the "average" 
pulsar is according to the "standard" spin-down law with n 5! This value may suggest the 
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importance of multipole components of pulsar magnetic field, or the deviation of the angle 
between pulsar rotational and dipole axes from n/2 (Manchester & Taylor 1977; Contopoulus 
& Spitkovsky 2005). From (8) we may easily get the relation between i) tr and v tr 

Vtriytr) = «C»(-i) tr ) 2 "», (9) 

which is shown in Figure 6 as a thick dashed line. The amplitude of the i)(t) oscillations, 
Aj>(t), may be easily computed from (7) using fits for i) + and i>_ in Figure 1, 

A i> (t) = Am 2 (t) = \\K(t)-i>-(t)} (10) 

The result is shown in Figure 6 as a solid line. The amplitude A„ is decreasing almost 
linearly with the increase of u tr , i.e. life time. This may be related to the decrease of either 
A{t) or Q(t). 

From Figure 6 it is seen that our simple model formulated in (2) - (7) describes the 
observational data rather well, at least for the pulsars with v > — 10~ n s~ 2 , i.e. with 
r ch > 10 4 years. For v < —7- 10~ 12 s~ 2 , the A„ is smaller than i) tr , and during the oscillations 
v > always, that explains the absence of pulsars of the negative branch of v — u diagram 
in this region. The situation changes dramatically for v > 10~ n s~ 2 , i.e. r ch > 10 4 years. 
Here A„ becomes larger than i) tr and pulsars with i) < and n < appear. In other words, 
the cyclic behaviour according to the model described by (2) - (7) begins to manifest itself 
in the observing quantities, - v and n begins to deviate significantly from the predictions of 
the "standard" spin-down law evolution. 

On the basis of these assumptions it is possible to get a rough estimate of the maximal 
period of the cyclic behaviour as the several times of the time scale of the vanishing of v 
for the young pulsars with the measured 'v as T rnax « —Av/'v. For PSR B1509-58 (with 
v « 2 • 10~ 21 s~ 3 and i> w -1.3 • 10~ 31 s~ 4 ) (Livingstone et al. 2005) T max ss 2 • 10 3 years 
(n min w 10~ 10 s" 1 ). For Crab (with v w 1.2 • 10~ 20 s~ 3 and v w -1.4 • 10~ 30 s~ 4 ) (Scott, 
Finger & Wilson 2003) T max w 10 3 years (VL min w 2 • 10~ 10 s" 1 ). On the other hand, the 
estimations of VL may be done using the relations (2) - (4), for the pulsars with several 
measurements of rotational parameters at different epochs with enough separation. So, the 
lower limit is Vt min ^/|Ai>/Ai/|. From the results of observations of the old (r oc 10 6 years) 
pulsars, B0823+26, B1706-16, B1749-28 and B2021+51 with 7 to 10 years separation, we 
have obtained the same estimate of T max m 10 3 -j- 2 • 10 3 years. Finally, the period of cyclic 
oscillations is likely to be in the range of 100 - 1000 years (2 ■ 10~ 10 < f2 < 2 • 10~ 9 s _1 ). 
Here we have suggested that the minimal cycle period has to be several times the pulsar 
observation period - 30 years. 
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Taking into account that A„ = Q 2 A is changing several orders of magnitude, our initial 
assumption of relative constancy of Q and significant decrease of A during the evolution is 
somehow justified. Assuming Q 10~ 9 s _1 , we get for the amplitude of frequency oscillations 
A = Ay/VL 2 values from 4 • 10~ 4 s _1 for the youngest pulsars till 1CT 10 s _1 for the eldest. 

4. Discussion and conclusions 

The physical reasons of the discussed cyclic variations of the pulsar rotational frequency 
may be similar to the ones of the well-known red timing noise. In the first place, the free 
precession (Stairs, Lyne & Shemar 2000) is the periodic effect most suitable for the proposed 
harmonic model. Note that the interaction of the superfluid core and the crust of the neutron 
star does not depress its development (Alpar 2005). Several processes had been proposed for 
the explanations of the red noise (Cordes & Helfand 1980; Cordes & Greenstein 1981; Cheng 
1987) - from the collective effects in the superfluid core till the electric current fluctuations 
in the pulsar magnetosphere. The possibility of its utilization for the explanation of the 
frequency variations on the time scale of hundreds of years is yet to be analyzed. 

The argument towards the similarity between the discussed variations and the red timing 
noise is the coincidence of the red noise v amplitude extrapolated according to its power 
spectrum slope (Baykal et al. 1999) to the time scale of hundreds of years, with Aj> in our 
model for the same z>, i.e. the same ages (see Fig. 6). 

Finally, note that the real evolution of the pulsar rotational frequencies may be sig- 
nificantly more complicated than the proposed simple model. The variations may be non- 
harmonic, their amplitudes, phases and frequencies may change randomly. The principal 
point is that all the pulsars evolve this way, and the time scale of such variations must sig- 
nificantly exceed several tens of years. That explains the anomalous values of the observed 
v and braking indices. 

This work has been supported by the Russian Foundation for Basic Research (grant No 
04-02-17555), Russian Academy of Sciences (program "Evolution of Stars and Galaxies"), 
and by the Russian Science Support Foundation. 
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Fig. 4. — The sample of pulsar cyclic behaviour on the i) — i> diagram. The solid line shows 
the evolution of a pulsar according to the simple harmonic model presented in (2) - (7) with 
the oscillation period T = 1000 years. The dashed line shows the corresponding secular 
evolution term i^ r (t), and the oscillations are around it, symmetric in linear scale. For 
v < —7 ■ 10~ 12 s~ 2 , the amplitude of the oscillations is smaller than the secular term, thus 
v is always positive, for larger values it changes its sign repeatedly, spending approximately 
half of the time in the negative region. 
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Fig. 6. — The i) — v diagram with the simple harmonic oscillations model. The solid line is 
the amplitude of the frequency second derivative variations according to (10), the dashed 
- the secular term i) tr from (9), and the dot-dashed lines are the envelopes of the oscillations 
i) tr ± Ay with l-o" ranges (dotted lines). Each pulsar spends the majority of its life time at 
or very near the envelopes. 



